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Abstract. - We study experimentally the dynamics of a dense intruder sinking under gravity inside 
a vibrated 22? granular packing. The surrounding flow patterns are characterized and the falling 
trajectories are interpreted in terms of an effectivive friction coefficient related to the intruder mean 
descent velocity (flow rules). At higher confining pressures i.e. close to jamming, a transition 
to intermittent dynamics is evidenced and displays anomalous "on-off" blockade statistics. A 
systematic analysis of the flow rules, obtained for different intruder sizes, either in the flowing 
regime or averaged over the flowing and blockade regimes, strongly suggest the existence of non- 
local properties for the vibrated packing rheology. 



Introduction. — During the last decades, several 
studies have shown amazing analogies among variety of 
very different systems pQ (colloids [2], foams [3], granu- 
lar materials [HE], supercooled liquids [E]). Indeed, close 
to the jamming transition these systems share common 
dynamical behaviors such as a glassy phase, aging, mem- 
ory or intermittency [6j. Based on these analogies, several 
propositions emerged aiming to describe and possibly un- 
derstand, in a general and unified way, this phenomenol- 
ogy and the associated rheology [7][S] . In the case of dense 
granular packing, the rheology under shear is described 
with good accuracy by local relations between a dynamical 
friction coefficient and a local dimensionless number, in- 
volving shear and confining pressure [4] : the inertial num- 
ber. However, the boundary between fluid-like and solid- 
like phases is still poorly understood. Experiments on 
avalanche flows suggest that the rheology could indeed be 
non-local near arrest |9ll0j . Conceptually, for the rheology 
of vibrated granular packing, the most elaborated point of 
view stems from kinetics theory [TT] . Whereas this simple 
vision is reliable to describe dilute systems, it is severely 
challenged for dense systems [15r[T7, 29 . It has been sug- 
gested in this case, that the notion of thermal temperature 
could be replaced by the more general concept of " effective 
temperature" describing the structural evolution between 
blocked configurations [12] : however, practically, it is quite 
difficult to apply this concept in order to obtain a clear 
constitutive picture. In experiments, a complication usu- 
ally comes from the way energy is input in the packing 



i.e. from the boundaries. It often leads to very inhomoge- 
neous bulk agitations, both in time and space. The probes 
themselves are hard to decouple from the agitation energy 
input [13]fl4] , which may alterate the reliability of the con- 
stitutive relations hence extracted [TJ] . This was somehow 
circonviened recently by Reis et al. [15] and Lechenault 
et al. [16], who studied 21? horizontal packing driven ho- 
mogeneously from the bottom; same for Keys et al. [17] 
using an hydrodynamic fluidization technique. In these 
experiments, which approach jamming by increasing the 
packing compaction, strongly heterogeneous and coopera- 
tive displacement modes were evidenced, very reminiscent 
to those present in undercooled liquids or glassy systems 

Recently, Caballero et al. [19] studied a sonofluidized 
31? granular packing and monitored the motion of various 
intruders inside. Even at weak driving acceleration, the re- 
sulting friction coefficients decrease significantly with ac- 
celeration. Eventually, the Coulomb threshold vanishes 
and one obtains a linear relation between force and veloc- 
ity. These relations are strongly dependent on the shape 
and the size of the intruders. Other means of activation 
were proposed using shear bands noise as actuators [201121] . 
In these last cases, empirical relations between force and 
velocity were extracted and all these works concluded on 
the non-local character of the constitutive relations. How- 
ever, in all these cases, the intruder was buried partially 
or completely in the bulk, thus hindering a direct visual- 
ization of granular motion in the intruder surrounding. 
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In the present work, we seek to extract the rheological 
properties of a vibrated vertical granular packing by mon- 
itoring the vertical motion of cylindrical intruders. This 
approach is very analogous to bead rheometers techniques 
currently used in fluid mechanics. Note also some recent 
work along those lines on intruder penetration in 21? or 
32? packing but in the absence of vibration 25,26 . This 
kind of experiments were also performed in other systems 
like colloids [2], glasses [24] or foams [27] . 

The paper is organized as follows. First, we present 
the experimental setup, a 22? experiment realization of a 
vibrated granular packing oriented in the direction of the 
gravity. Then, the descent dynamics of heavy intruders are 
monitored in association with the resulting counterflow of 
the surrounding particles. We show that the analysis of 
the intruder trajectories may probe a rheological behav- 
ior via a relation between an effective friction coefficient 
and the mean intruder descent velocity. We analyze the 
intermittent statistics evidenced at higher depths as an 
alternation of flow and blockade. We finally question the 
existence of a local relationship accounting for the rheol- 
ogy of the dense granular medium under vibration. 

Set-up and experimental protocol. Figure [TJa) 
displays a snapshot of the experimental set-up. The cell 
is a vertical rectangular container of lateral size L — 275 
mm and thickness 3.5 mm slightly larger than the thick- 
ness of the disks forming the model granular packing. We 
use a bi-disperse mixture of around 1100 big (diameter 
df, = 5 mm) and 1600 small (diameter d s — 4 mm) ny- 
lon disks of thickness 3 mm. The intruder is a metallic 
cylinder of adjustable radius R and same thickness as the 
grains. The average intruder density can also be adjusted 
by drilling a hole at its center. The density contrast with 
the surrounding vibrated granular medium of density p is 
defined as ^ and is always chosen positive. Therefore, 
the intruder will sink due to buoyancy if the surrounding 
medium is fluidized enough. 

The way vibration energy is input in the system, dif- 
fers from what has been done in previous experiments. 
Here, we seek to avoid synchronization between the global 
motion of the packing and the sinusoidal motion of the 
driving plate. It is known that for accelerations larger 
and around gravity <?, the packing trajectory and the en- 
ergy input due to repetitive collisions with the bottom 
plate, are strongly synchronized, unless a very high accel- 
eration level is reached (up to 30g! 28 ). Note also that 
for driving accelerations larger than g, surface instabili- 
ties such as Faraday waves may also occur which render 
the packing spatially inhomogeneous [3D] , Here, the idea 
is to break as much as possible this spatio-temporal syn- 
chronisation responsible for the alternation of phases at 
drastically different agitation levels. To this purpose, we 
place a so-called " vibrating chain" at the bottom of the cell 
(see Figure Ufa) bottom) which consists of two arrays of 
large and dense disks activated from below by 24 electro- 
magnetic pistons oscillating vertically with an amplitude 
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Fig. 1: (a) View of the cell containing the bi-disperse pack- 
ing and a cylindrical intruder of radius R = 10 mm with 
p = 0.62 g.cm~ 3 '. The dashed rectangle is the visualization 
field, (b) Vertical positions of the intruder as a function of 
time for 10 experiments with identical conditions of vibra- 
tions (radius R = 10 mm, density contrast — 2 = 16.15). 
Inset : vertical increment of displacement for a time lag Is 
for one of the intruder trajectories presented in the main 
graph (same color). 

A p inside the granular packing and thus, transferring vi- 
bration energy to this " chain" . A detailed study of bulk 
agitation and packing density hence obtained at steady- 
state, is discussed elsewhere [31] , 

In this report, the packing is vibrated under fixed driv- 
ing conditions : piston amplitude Ap — 6 mm and driving 
frequency fjj = 20 Hz with a phase shift of 7r between 
neighboring pistons. We measured for this system a driv- 
ing RMS acceleration of 2g in the vertical direction. In a 
typical experiment, the vibration is activated 15 minutes 
before we place the intruder at the piling surface. Its de- 
scent is monitored by a CCD camera taking pictures every 
second. For each intruder we perform 10 experiments in 
similar experimental conditions. 

Intruder descent dynamics. — Figure [ljb) shows 
the vertical positions z(t) of an intruder of radius R = 10 
mm as a function of time. The first noticeable feature is 
that the intruder motion is quite irregular with alternate 
phases of flow and blockade. The deeper the intruder is in 
the packing, the more likely blockade seems to occur (see 
inset of figure [ljb)). This intermittent dynamics yields a 
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Fig. 2: (a) Average displacement field for a time lag corresponding to an intruder displacement over its radius 
R = 10 mm and for a coarse-graining size AW = d s . The color code represents the amplitude of the coarse-grained 
displacement normalized by the intruder radius and the vectors indicate the direction of grain motion. Note the 
horizontal distance l v between the intruder center and a vortex center, (b) Horizontal positions of left and right vortex 
centers as a function of time ( density contrast —^=16.15), ((•) right vortex and (•) left vortex), (c) Mean horizontal 
distances l v between the intruder and the vortex centers as a function of the radius R (symbol colors as in (b)). 



large distribution of arrival times at the bottom of the cell 
as observed from the different vertical trajectories taken 
in identical conditions. 

Figure [^a) displays the coarse-grained displacement 
field over a time lag r defined using the mean intruder 
sinking time over a distance R. A feature of these displace- 
ment fields is the existence of two vortices accompanying 
the intruder's descent and corresponding to recirculations 
of grains around the intruder. Note that during a descent, 
one of the vortices seems to be dominant in amplitude but 
it can be located on either side of the intruder. In figure 
[2fb), the distances l v joining the intruder's center to both 
vortex centers are represented as a function of time. From 
Figure[5fb) , we see that during the whole fall both vortices 
seemed to stay at a fixed distance l v (R) from the intruder. 
We thus represented this mean distance l v (R) for experi- 
ments done with different intruder sizes and observed that 
these distances scale almost linearly with the radius R ( 
Fig. ^c)). 

Effective friction. — One objective of this study is to 
shed light on the rheological properties of granular mat- 
ter under vibration. In either static, quasi-static or dense 
flow situations, the central constitutive parameter describ- 
ing the granular rheology is the ratio p between shear 
stress and confining pressure (see 0] and refs. inside), 
i.e., a Coulomb friction coefficient. Scaling analysis shows 
that for an intruder of radius R, the driving force per 
unit of thickness Fq is set by gravity, thus Fq cc ApgR 2 . 
This force must be balanced at steady state by the re- 
sulting friction force exerted on the intruder boundary 
Fs oc ppgzR. Thus, a relation can be obtained for the 
friction coefficient p which scales like p oc ApR/pz, where 
—P- is the density contrast, R the intruder radius and z its 
vertical depth. This argument can be made more quantita- 
tive in the framework of a simple reference model involving 
a constant dynamical friction p so that effective friction 



and dynamical friction would match (p e = p). The force 
balance can be computed exactly by considering a local 
confining pressure imposed by gravity P — pgz and by 
integrating the friction force around the intruder. This is 
the base for our definition of an effective friction coefficient 
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This relation motivates a parametric study of the in- 
truder motion by varying its size and its apparent density. 
For each realization i, we measure the time Ti(z) spent 
by the intruder to span the interval [z — Az/2, z + Az/2]. 
Then, the mean descent velocity is defined as V(z) =< 
Az/Ti(z) >i, the average being taken over 10 realizations. 

Dimensionless number. — In figure^a), we display 
for an intruder of radius R = 10 mm, the mean descent 
velocity V(z) obtained for Az = 20 mm and for different 
density contrasts. As expected the mean velocity increases 
with density contrast. We also note the very regular de- 
crease of mean velocity with depth z. In order to rescale 
the different curves, we define for each depth (z), a lo- 
cal mean Froude number as Fr — V / \fP~Jp = V/y/gz. 
In Figure EJb), we observe a rescaling of all the data on a 
common flow rule, which justifies the existence of an effec- 
tive friction coefficient describing the rheology of the gran- 
ular intrusion. Interestingly, this flow rule is of the hard- 
ening type, i.e., the effective friction increases with the 
local drag velocity, a behavior also observed for sheared 
dense granular media [4 without [26] or with vibration 

HSH2U. 

Blockade and flow statistics. We previously no- 
ticed (see Figure [Tib)) that a prominent feature of the 
intruder descent dynamics is the intermittent sequel of 
blockades and flows. This clearly has a strong impact on 
the effective mean velocity as measured previously. In the 
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Fig. 3: (a) Mean velocity V(z) as a function of depth 
z ± Az/2 with Az = 20 mm and for R = 10 mm and 
different density contrasts (b) Effective friction /^(z) 

— for different 



as a function of Froude number Fr(z) 
density contrasts. 



following, we will seek to separate and characterize both 
regimes independently. The method we design to extract 
the blocked phases is based on the use of a displacement 
threshold Sz at a time scale 6t = 1 s below which we con- 
sider the intruder is blocked (see the dashed line in inset of 
figure [2(b)). This choice is somehow arbitrary; however, 
we rationalized it by defining, for each descent curve, a 
value such that the jump distribution in the blocked phase 
is close to a symmetric Gaussian curve with a zero average. 
For each trajectory, we identify the depths at which the 
intruder gets blocked and starts flowing again as well as 
the time spent in both phases respectively. Figure 0] dis- 
plays the distribution of times spent in the blocked phase 
tb at different depth values (z ± Az/2). The striking fea- 
ture is the outcome of a power-law distribution t~ a , with 
an exponent smaller than 2 (a = 1.5). Such a distribution 
is the sign of an anomalous statistics as no typical time 
scale can be defined. To characterize the sojourn in the 
blocked phase, we thus consider the fraction of time $6 
spent in the blocked phase at different depths. 

In figure El is plotted as a function of the effective 
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Fig. 4: Probability density function of blockade times 
n(t(,) for R = 10 mm and — 16.15, at different depths. 




Fig. 5: Fraction of time $b(z) spent in the blocked phase 
as a function of the effective friction \x & (z) for R = 10 mm 
and different density contrasts — . 



friction ji e {z). At low $(, is large and almost constant. 
At fi e = 0.7 a sharp decrease of the fraction of blocked time 
is observed. This transition from blockade to flow suggests 
a generalization of the static Coulomb threshold in the 
presence of vibration. Note that this effect of threshold 
distribution was also evidenced in numerical simulations 
of sheared granular packings under random noise |29) . 
Now we want to characterize the flowing phases. We de- 
fine a flowing time variable £f by removing sequentially 
the blockade periods of the trajectories. Figure EJa) dis- 
plays the flowing trajectories i.e. the intruder positions as 
a function of tp. A clear collapse of all the trajectories is 
observed and an average is performed over 10 experiments 
done in identical conditions in order to define a mean flow- 
ing trajectory z(tp). A flowing velocity is then computed 
at each depth Vp — This procedure allows to define 

a new Froude number Frp = j£= , (with the subscript 



F for flow) that we call the dynamical Froude number as 
it applies to the flowing phase only. 
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In figure [U^b), we display the relation fj, e (Frp) for an 
intruder of radius R — 10 mm and for various density con- 
trasts. Interestingly, we also observe a data collapse onto 
a unique curve. 
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Fig. 6: (a) Vertical positions z of the intruder with R = 10 
mm corresponding to figjljb), as a function of the flowing 
time variable tp. (b) Dynamical flow rule : effective fric- 
tion fJ, e (z) as a function of the dynamical Froude number 

FVjp{z^ — f' 11 ' \':>vwiiw ili'iwiiv (villi < ^LL 



for various density contrasts =£. 



Towards a non-local rheology ?. — To extract a 
local rheology that could reflect constitutive properties 
of granular matter under vibration, we must at the end, 
provide a relation between a constitutive coefficient (fric- 
tion) and local macroscopic fields such as pressure or shear 
rate and this, independently of the intruder's presence. 
To address explicitly this question, we computed system- 
atically the mean and dynamical flow rules correspond- 
ing to intruders whose diameter varies between 3.75 and 
6.5 small granular sizes d s . We obtained in each case - 
like what was presented previously for R = 10mm - a 
rescaling of the data with different density contrasts. By 
visualization of the intruder flows (see Figures EJa) and 
(b)), we have seen that shear is exerted on the intruder 
side at a distance of the order of l v (R) yielding a typi- 
cal shear rate 7 = TTE)' Following the analysis provided 
for rapid dense flows, we may try to transform the flow 



rules into a local rheology law using an inertial parameter 
j = _jd^ = dsVizi = d^p If j j rheol applies, 

a relation fi e (I) should show up independently of the in- 
truder size. We see in Figures [Tfa) and (b) that the flow 
rules for different intruder's size can be rescaled by divid- 
ing the Froude numbers by a value Fr*(R) for the mean 
flow rule or Fr* F (R) for the dynamical flow rule. These 
rescaling values are presented as a function of the intruder 
radius R in insets of Figures[7Ja) and (b). However, we ob- 
served that for both flow rules the corresponding lengths 
Ft (R)d s and Fr F (R)d s , grow much faster with R than 
l v (R) . Therefore these rescalings could not provide a value 
such that the Froude ratios could be identified with an in- 
ertial parameter independent of R. Thus, it becomes clear 
that contrarily to the situation of rapid dense flows, the 
flow rule cannot be transformed into a local rheology us- 
ing an inertial parameter. Whether this is the sign of true 
non-locality in the constitutive relations associated with a 
growing length scale near jamming or the unravelled in- 
fluence of a supplementary field such as "fluidity" or " 
granular temperature" is left for future discussions. Note 
finally that under constant velocity and confining pres- 
sure, the effective friction increases when the intruder size 
decreases. This is the so-called "geometrical hardening" 
effect already noticed by Caballero et al. |19j . 

Conclusion. — In this experimental study, we mea- 
sured the dynamics of dense intruders sinking in a vibrated 
granular packing. The flow fields around the intruder was 
measured and characterized by two side vortices accompa- 
nying the intruder's descent. The mean position of these 
vortices scales with the intruder size. The sinking dy- 
namics is interpreted as as an effective friction coefficient 
increasing with the mean falling velocity. For a given in- 
truder size, when the mean falling velocity V(z) is rescaled 
as a Froude number Fr(z) — , where P(z) is the 

\/P{z)/P 

local confining pressure, a rescaling of all the flow curves 
is evidenced, for different intruder densities. This defines 
the so called " mean flow rule" characterizing the intruder 
rheology. An important feature is that at large confining 
pressures, the intruder dynamics is intermittent and dis- 
plays alternating phases of blockade and flow. The time 
spent in the blocked phase follows a statistics with a decay- 
ing power-law displaying an anomalous exponent around 
1.5. This is quite reminiscent of the "on-off" intermit- 
tency process [32] where the blocked phases show up with 
a similar residence time distribution, also evidenced for 
the flow out of a vibrated hopper [33]. We interpret the 
transition between blockade and flow as an effective fric- 
tion threshold appearing as a " noisy" generalization of the 
Coulomb threshold. When isolated, the dynamics of the 
flowing phases, displays a similar data collapse defining a 
" dynamic flow rule" , however with a different functional 
form as for the "mean flow rule". Finally, for different 
intruder sizes, we could collapse both flow rules using a 
rescaled value of the Froude number increasing exponen- 
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Fig. 7: (a) Mean flow rules for different intruder sizes rep- 
resented as a function of a rescaled Froude number -Sr 

Ft 

for different intruder radius R (fit A$.x Al with Aq = 0.56 
and Ai = 0.29). Inset : rescaling Froude number Fr as 
a function of the intruder radius, (b) Dynamical flow rule 
obtained in the flowing phase as a function of a rescaled 
Froude number ^4- (fit A .x Al with Aq = 0.56 and 

A\ = 0.37). Inset : rescaling Froude number Fr F as a 
function of the intruder radius. Note the log scales for the 
vertical axis. 



tially with the intruder size. We discuss the fact that this 
very strong dependance might be an indication that the 
effective rheology of the agitated packing is non-local. 
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